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Chip-based photonic quantum computing is an emerging technology that promises much speedup over 
conventional computers at small integration volumes. Particular interest is thereby given to 
polarisation-encoded photonic qubits, and many protocols have been developed for this encoding. However, 
arbitrary wave plate operation on chip are not available so far, preventing from the implementation of 
integrated universal quantum computing algorithms. In our work we close this gap and present Hadamard, 
Pauli-X, and rotation gates of high fidelity for photonic polarisation qubits on chip by employing a 
reorientation of the optical axis of birefringent waveguides. The optical axis of the birefringent waveguide is 
rotated due to the impact of an artificial stress field created by an additional modification close to the 
waveguide. By adjusting this length of the defect along the waveguide, the retardation between ordinary and 
extraordinary field components is precisely tunable including half- wave plate and quarter- wave plate 
operations. Our approach demonstrates the full range control of orientation and strength of the induced 
birefringence and thus allows arbitrary wave plate operations without affecting the degree of polarisation or 
introducing additional losses to the waveguides. The implemented gates are tested with classical and 
quantum light. 



Quantum information science addresses various fundamental questions on how to harness quantum 
mechanical phenomena for storing, processing and transmitting information 1 . Therefore, quantum 
information is both fundamental science and a progenitor for new technologies 2 . One of the most 
intricate future technology is the quantum computer which promises exponentially faster operation for particular 
tasks such as factorising, searching databases and simulating important quantum systems 1 . Along these lines a 
number of physical systems for the development of well -performing quantum computers have been investigated, 
such as ultracold atoms and trapped ions, spin chains, cavity quantum electro dynamics, superconductors, and 
liquid state NMR 3 " 8 . In addition, those systems involving photonic quantum states seem likely to play a central 
role as light is a logical choice when dealing with quantum communication, metrology and lithography. 
Moreover, photons are in general more immune to decoherence than other quanta. 

However, a downside of photonic quantum architectures is the complex setting, at least when using bulk 
optical systems: Building more complex networks of logic gates encounters severe size and stability limitations. A 
promising approach to miniaturising and scaling optical quantum circuits is to use on-chip integrated photonic 
waveguides 2 , which show strong improvements in performance due to high stability and low noise - crucial 
aspects for high-fidelity classical and quantum interference. 

In recent years various breakthroughs were presented, in which the integration of quantum circuits using 
optical waveguides played an important role. The first landmark was the implementation of a photonic quantum 
gate on a chip, employing silica-on-silicon waveguides 9 . Using this approach it was shown that it is possible to 
realise one of the most important quantum gates - the C-NOT gate - on a chip, for path-entangled photons. This 
first seminal work was followed by various follow-up publications, in which the potential of integrated quantum 
optics was explored: Integration of controlled phase shifters in integrated interferometers has been used to 
manipulate entangled states of up to four photons and to demonstrate on-chip quantum metrology 10 . A compiled 
version of Shor's quantum factoring algorithm was performed on chip consisting of several one- and two-qubit 
gates 11 . Moreover, an on-chip waveguide architecture was employed for the demonstration of continuous-time 
quantum random walks of path-entangled photons 12 . And only recently even a boson sampling computation was 
executed (besides using bulk optics 13 ) on a photonic chip 14 " 16 . 
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Importantly, photons exhibit polarisation, a two -level degree of 
freedom that carries numerous benefits when encoding qubits 17 . 
Using polarisation entanglement is highly beneficial as various 
entangled-photon sources are available 18 , and many quantum 
information protocols and experiments have been developed for this 
encoding 19 22 . An arbitrary polarisation transformation, which 
means to have the ability of preparing a single qubit in any desired 
state, was shown for telecom wavelength using the electro -optical 
effect in titanium- diffused lithium niobate waveguides 23 ' 24 . Pushing 
towards the visible regime it turned out that waveguide architectures 
fabricated using the femtosecond laser- writing technique 25 are par- 
ticularly useful for polarisation-encoded on-chip quantum optics 
due to the low birefringence and, consequently, low decoherence 26 . 
This sophisticated fabrication approach allowed the realisation of C- 
NOT gates with polarisation encoding 27 . Laser-written waveguides 
are also an exceptional platform for implementing discrete-time 28 
and continuous -time quantum random walks 29 . Moreover, three- 
dimensional geometries are feasible with such architectures 30 . 

In our work, we go the next consequent step in the development of 
laser-written integrated quantum circuits and present photonic 
Hadamard gates, Pauli-X gates, and rotation gates on chip. In par- 
ticular the Hadamard gate is one of the three universal quantum gates 
(besides the C-NOT gate and the 71/8-gate) which, according to the 
quantum circuit model, represent a set of building blocks that effi- 
ciently approximate any conceivable quantum operation if arranged 
in a given finite sequence 1 . Therefore, our results constitute a crucial 
leap towards polarisation- entangled quantum computing on chip. 

Results 

For the fabrication of the integrated quantum gates as well as the 
general wave plates we employ the femtosecond laser-writing tech- 
nique 25 . Only recently it was shown that the birefringence in a laser 
written waveguide can be precisely tuned by employing particular 
line defects close to the waveguide 31 . However, in this work the 
optical axis of the birefringent waveguides was always fixed either 
along the vertical or the horizontal direction, preventing an actual 
Hadamard operation. We generalise this approach in order to freely 
choose the orientation of the waveguide's optical axis, whereas the 
waveguides themselves remain in their original orientation. With our 
approach we are able to implement arbitrary desired polarisation 
transformations, including the Hadamard operation. 

Importantly, laser-written waveguides are surrounded by an 
extended stress field that is created during the writing process 32 . So 
far, this stress field has been used to alter the strength of the birefrin- 
gence of adjacent waveguides 31 ; however, there was no impact on its 
vectorial character. 

Our experimental analysis of this phenomenon delivers the sur- 
prising result that, when an additional waveguide is prone to such a 
stress field, an artificial birefringence is induced in this guide, which 
results in a reorientation of the optical axis as a function of the 
relative position of the two guides. As the relative position (in dis- 
tance and angle) can be tuned with high precision, strong tilting of 
the birefringent axis can be implemented at will, allowing arbitrary 
desired wave plate operations on the light that propagates in the 
waveguide exposed to the stress field. A sketch of this effect is shown 
in Fig. 1(a). To prevent the system from coupling light to the defect, 
the defect is written slightly above the destruction threshold. The 
main advantage of our approach is that the original waveguide shape 
remains unchanged irrespective of the optical axis' orientation that 
causes the wave plate operation. Hence, there is no additional friction 
of the propagating light caused by the artificial birefringence The 
additional losses in presence of the defect are <0.01 dB/cm. 

In order to test the functionality of our approach, we first perform 
a series of classical measurements using free space cw laser light at 
X = 815 nm (Spectra Physics Tsunami) that is launched into the 
integrated wave plate using a 10 X objective (NA = 0.25) and 




waveguide 



Figure 1 | Settings, (a) Idealised schematic showing the cross section of the 
waveguide arrangement where additional stress fields induce a 
reorientation of the waveguide's optical axis, (b) Sketch of the writing 
setting, with which the quantum gates are fabricated (here for 3 different 
orientations of the defect relatively to the waveguide). Higher pulse 
energies and lower writing velocity for the defects causes their missing light 
guiding properties and, hence, prevent loss due to coupling. 

projected onto a power meter using a 4X objective (NA = 0.10). 
In our experimental analysis the distance between the waveguide and 
the defect is r = 20 urn, such that the introduced stress field is large 
enough to change the optical axis of the target waveguide without 
disturbing its mode profile. As a benefit it turned out that only a 
defect on one side is necessary to obtain the desired effect while the 
degree of polarisation is still retained (>1000 : 1). We characterise 
the strength of the birefringence as well as the orientation of the 
optical axis following a method proposed by Yang et al. 33 . Two 
crossed polarisers encircle the optical chip and are rotated simulta- 
neously over the desired range. The signal on the power meter drops 
to zero only if the polarisers are orientated along the fast or slow axis 
of the investigated wave plate. Further rotating the polarisers by 45° 
the maximum power is transmitted and discloses the strength of the 
birefringence. Detailed information can be found in the method 
section. 

We extract a birefringence in the waveguide that is An = 0.65 ... 
1.32 X 10" 5 , depending of the azimuthal angle of the defect caused by 
the resulting overlap of the stress fields. The orientation of the optical 
axis as a monotonous function of the azimuthal angle 0 of the defect 
is plotted in Fig. 2(a). For increasing 0 also the tilt of the waveguide's 
optical axis a increases (in a nonlinear fashion) and reaches a 
maximum of a max = 90° for 0 = 90°. In order to understand our 
experimental data, we developed a simple heuristic model. As the 
birefringence in the waveguide is increased by approximately one 
order of magnitude in the presence of a defect, we consider only the 
defect's stress fields, represented by their centroids on either sides of 
the defect with coordinates {r b 0 t ) with respect to the waveguide. 
Assuming that the impact s f of such a stress field i on the optical axis 
decreases exponentially with growing distances by the decay rate p, 
we obtain for our data the fit function 

■W-TTT-- ' (1) 

i = l,2 f = l,2 



with Yi = y (r sin 0 + d yeYt ) 2 + (r cos 0 + ( — l) J dhor) and dhor.ven as 
horizontal or vertical deposition, respectively, of the stress fields' 
centroids compared to the centre of the defect. The best fit shows 
p = 20(3), um d hor = 7.95(10) um, and d vert = —4.2(6) um, sug- 
gesting that the centroids tend to lie below the envisaged level, which 
is caused by an elongated shape of the defect leading to an offset of 
0 O = 11.9° for zero- crossing. For smaller angles the experimental 
data deviate from the predicted line leading to a tilt of the wave- 
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Figure 2 | Measured wave plate operations, (a) Experimental data and best fit model of the reorientation of the optical axis as a function of the 
azimuthal position of the defect. Larger error bars result from near X operation of the integrated wave plate, (b) Intensity ratio P ± /P to tai between the power 
of the perpendicularly polarised light (H or V polarisation) and the total power as a function of the interaction length between defect and waveguide 
(assumed to be proportional to A(p). The angle of the optical axis is fixed to a = 45° (black dots) or a = 22.5° (red triangles) and P±/P to tsi should follow 
(sin(2a) sin(A^/2)) 2 . 



guide's optical axis bya = — 6° at 0 = 0°. This is due to the sequenced 
writing of the waveguide followed by the defect that has to be written 
in or through a yet present glass modification, which affects most 
notably the stress field on the left side of the defect. Hence, both stress 
fields will lie in different depths and the assumed symmetry, d ver t,i — 
^vert,2> is broken. 

Clearly, with this tuning degree of freedom we can produce arbit- 
rary wave plate operations, including the Hadamard gate, the Pauli-X 
gate, and a rotation gate. First, the realisation of those gates depends 
on the correct orientation of the optical axis (a = 22.5° or 45.0°) 
which is tuned by the position of the defect given by the angle 0. 
Second, it depends on the defect length z, since for increasing inter- 
action length the phase delay Ag? between the slow and the fast axis 
increases as well (see Fig. 2(b)). Table I outlines the corresponding 
parameters for implementing the gates. 

The fabricated Pauli-X and Hadamard gates are particularly inter- 
esting for quantum applications, as especially the latter is part of a 
simple set of universal quantum gates (besides the C-NOT gate and 
the 7i/8-gate). Hence, we experimentally demonstrate their function- 
alities also in the single photon regime. When a single photon in 
either horizontal \H) or vertical | V) -polarisation representing the 
1 0) or 1 1) state, respectively, enters a single qubit gate, the gate opera- 
tion is generally described by the basis transformation 



\V)=a\ti) + b\V)=(l)^\V oat )=M\V lD ) 



(2) 



where M is a 2 X 2 matrix. In case of the Hadamard gate, the 
transformation matrix 



M Ha d = 



Mi 



(3) 



is applied, which transforms a photon in the \H) state into the (diag- 
onal) \D) state 



Table I | Gate parameters 



gate 


a 


0 


A(p 


Z 


Hadamard 


22.5° 


20.5(2)° 


71 


3.6(1) cm 


Pauli-X 


45.0° 


29.7(3)° 


71 


5.6(1) cm 


rotation 


45.0° 


29.7(3)° 


% 


2.8(1) cm 



Outline of the corresponding parameters for implementing the Hadamard, Pauli-X, and a rotation 
gate. Note, only the relative positioning angle of the defect 9 and the defect length zare 
experimentally tuned. 



\H) 



\D) = 



\H) + \V) 
V2 



and a photon in the | V) state into the (anti-diagonal) | A) state 

|H>-|V> 



|V>-|A> = - 



V2 



(4) 



(5) 



Discussion 

The performance of the fabricated Hadamard gate is best reflected by 
its measured matrix that is shown in Fig. 3(a) and which depicts the 
response of the system as a function of the input state. It is clearly 
seen that indeed \H) and | V) input states are successfully transformed 
into |D) and |A) states, respectively. From our measurement data, we 
retrieve the experimental transformation matrix in the \H), | V) basis 
for our Hadamard gate 

1 / 1.014-/0.012 0.996 

V^V 0.986 -1.003-/0.025. 



(6) 



which is very close to the ideal (theoretical) Hadamard operation 

withafidelity (^*> = |(»P|M* ieal M exp |T)| 2 )'of^ |1 ' > =0.999(5). 

The error of the matrix elements ranges from 0.010 to 0.016 and 
originates from the Poissonian statistics of the single photons and 
the precision of the adjusted measurement bases. From the modulus 
of the matrix elements we deduced the orientation of the optical axis 
to be a Ha d = 22.3(4)°. The optimal value of a = 22.5° describes a 
transformation, e.g., from \H) to \D), i.e., with equal coefficients for 
\H) and | V). The slight discrepancy suggests that |H) and | V) are not 
perfectly balanced. The non-zero imaginary part indicates deviations 
from a perfect half-wave plate operation, where A(p should reach 
exactly n (see Table I). Hence, depending on the portion of the 
imaginary part the resulting state suffers from left or right rotating 
components. The experimental gate shows an almost vanishing ima- 
ginary part resulting in an A^ Ha d = 1.01(2)71 fulfilling the desired 
condition. 

The Pauli-X gate is defined by the basis transformation 



M PX = 



i.e., it flips the |H) and \ V) states: |H) |V) and vice versa. The 
experimental Pauli-X matrix (see Fig. 3(b)) reads as follows 
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Figure 3 | Measured gate matrices. Real and imaginary part of the 
measured transformation matrices for the Hadamard operation (a) and 
the Pauli-X operation (b). 



M exp - 

iVi px — 



-0.017-/0.096 0.997 
0.995 0.031-/0.070 



(8) 



and yields a fidelity of J 7 ^^ = 0.992 (7) . The extracted values of the 
wave plate operation are a PX — 45.7(7)° and A^ PX = 1.05(2)71 , which 
confirms the first implementation of a polarisation swapping device 
on a chip. 

In conclusion, we demonstrate a new approach to fabricate arbit- 
rary wave plates in optical integrated devices using induced birefrin- 
gence allowing precise control of the waveguide's optical axis 
orientation. Therewith, we successfully performed a Hadamard 
and a Pauli-X operation on both classical and quantum states of light, 
which constitutes a benchmark for a full manipulation of polarisa- 
tion-encoded qubits on-chip. Our results disclose new scenarios for 
integrated quantum computation, quantum information processing, 
and quantum simulation, and possibly find application in standard 
photonic devices whenever integrated and finely controllable 
birefringent devices are of interest. 

Methods 

Direct waveguide writing. We write our waveguides into transparent fused silica 
wafers (Corning 7980 ArF Grade), using ultrashort laser pulses (t ~ 150 fs, 
X = 800 nm) that are focused 250 am below the sample surface using a 20 X 
microscope objective (NA = 0.35). The defect that rotates the optical axis of the 
original waveguides is written with a higher pulse energy and lower writing velocity 
slightly above the destruction threshold, such that these defects do not guide light. 
However, their stress field induces a strong birefringence in the waveguide, with an 
orientation that is determined by the relative position of the defect with respect to the 
guide. For a sketch of the writing setting see Fig. 1(b); a microscope image of a written 
structure is shown in the inset of Fig. 2(a). The actual writing speed is 60 mm/min at a 
pulse energy of 200 nj and a repetition rate of 100 kHz (Coherent Mira/Reg A). Such 
waveguides exhibit low propagation losses (<0.3 dB/cm) and almost vanishing 
birefringence (Ahh,v ~ 10~ 6 ). With a supported mode field diameter of 12 urn X 
15 um coupling losses of 3 dB are obtained with standard single mode fibres. 

Determining orientation and strength of the birefringence. The birefringent 
properties of the individual waveguides are obtained by placing a polariser (Pol) at 
each side of the sample. The polariser after the quarter- wave plate (QWP) is used for 
preselecting the polarisation direction of the light beam at the input facet. The crossed 
polariser at the output is then used to analyse the transmitted light. From the positions 
of the transmission minima as a function of the polariser's rotation angle (see 
Fig. 4(a)) and the contrast between maxima and minima we are able to deduce the 




MM fibre 



glass sample 



HWP 2 



Figure 4 | Characterisation setups, (a) Classical characterisation of an 
arbitrary wave plate operation. The relative output of crossed polarisers 
gives information on the birefringence of the embedded waveguide, (b) 
Setting for the characterisation of the integrated quantum gates using 
single photons. One photon of a pair heralds the other one, which is in this 
case feed to the corresponding quantum gate via a butt-coupled fibre. The 
HWP 2 determines the measurement basis to be I H) and I V) or ID) and 
IA>. 

orientation of the optical axis in the waveguide as well as the strength of birefringence. 
The underlying formula reads as 



-^transm 
Jtotal 



[2(a Po i-a wp )] [1- cos(Ankz)} 



(9) 



with Jtransm/ftotai as relative intensity transmission, a Po i as orientation of the first 
polariser and oc wp as optical axis orientation of the evaluated wave plate, An as strength 
of the wave plate's birefringence, k as wave number and z as defect length. It can be 
seen that the transmitted intensity drops to zero for a Po i = a wp + m-90° (meZ). 
Finding those points reveals the orientation of the wave plate's optical axis. In 
contrast, the maximum intensity transmission always occurs at 45° away from the 
minimum position. With the used wavelength and the thickness of the wave plate the 
strength of birefringence can be extracted. 

Quantum light measurement. We generate single photons at X = 815 nm using a 
standard type I spontaneous parametric down- conversion source. A BiBO crystal is 
pumped by a 100 mW, 407.5 nm laser diode producing |H) polarised single photon 
pairs collected by polarisation maintaining (PM) fibres. Photon 1 is used to trigger the 
counts of the photon 2. This prevents the measurement from accidental counts 
caused by external photons or dark counts of the utilised avalanche photodiodes 
(APDs). By rotating half- wave plate (HWP) 1 photon 2 can be set to \H) or | V) before 
it passes the glass sample. The measurement basis can be aligned to |H) and | V) or |D> 
and \A) by rotating HWP 2 to 0° or 22.5°. Since a single photon state can be 
considered as a pure state, which means it is located on the surface of the Bloch sphere, 
the projections determined by HWP 2 yield in general two possible points | X F) A and 
| X F) B for the measured state. These points only differ in the complex phase distributed 
symmetrically around 0 or n, respectively. In other words | V F) A and | V F) B describe the 
same state irrespective if it is left or right circular. To get a knowledge about the sign of 
the complex phase and to have the ability to associate the output state to either | V F) A or 
|*F) B , a qualitative examination using a quarter- wave plate is necessary, which was 
already done by the classical experiment. 

The resulting state for a |H) and | V) input state is measured by collecting the 
outputs of a polarising beam splitter (PBS) for different bases via multi mode (MM) 
fibres. A sketch of the setting for the single photon characterization of the gates is 
shown in Fig. 4(b). 
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